Analysis of negative magnetoresistance. 
Statistics of closed paths. I. Theory 
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Statistics of closed paths in two-dimensional (2D) systems, which just determines the interfer- 
ence quantum correction to conductivity and anomalous magnetoconductance, has been studied by 
computer simulation of a particle motion over the plane with randomly distributed scatterers. Both 
ballistic and diffusion regimes have been considered. The results of simulation have been analyzed 
in the framework of diffusion approximation. They are used for calculation of the magnetic field 
dependence of magnetoconductance in the model 2D system. It is shown that the anomalous magne- 
toconductance can be in principle described by the well known expression, obtained in the diffusion 
approximation, but with the prefactor less than unity and phase breaking length which differs from 
true value. 

PACS numbers: 73.20.Fz, 72.20.Dp, 72.10.-d 



I. INTRODUCTION 

It is well known that the interference of electron waves 
scattered along closed trajectories in opposite directions 
produces a quantum correction to the conductivity. An 
external magnetic field applied perpendicular to the two- 
dimensional (2D) layer destroys the interference and sup- 
presses the quantum correction. This results in anoma- 
lous negative magnetoresistance, which is experimentally 
observed in many 2D systems. This phenomenon can be 
described in the framework of quasiclassical approxima- 
tion which is justified under the condition kpl ^> 1, where 
kp is the Fermi wave vector, I is the mean free path. In 
this case the conductivity correction is usually expressed 
through the classical quasiprobability for an electron to 
return to the aceaof the order Xpl (Xf — 2n/kp) around 
the start pointEfcl 
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where do = e 2 kpl/ (2ttK), and W stands for the 
quasiprobability density of return (quasi- means that W 
includes not only the classical probability density, but 
the effects of interference destruction due to an external 
magnetic field and inelastic scattering processes). In or- 
der to calculate a magnetic field dependence of negative 
magnetoresistance the quasiprobability W is represented 
as a sum of contributions of closed paths with N colli- 
sions, Wn- Then, Eq. (fy) can be rewritten as 
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where Go = e 2 / (2it 2 Ti). Here, only paths with N > 3 
are taken into account, because the paths with N = 1, 2 



have zero areas and their contributions are not influenced 
by the magnetic field. 

The expression (^J) for backscattering quantum correc- 
tion is true for an arbitrary magnetic field, any anisotropy 
of scattering and distribution of scatterers, and various 
relationship between phase and momentum relaxation 
times, t v and r, respectively. The sum (^) is usually cal- 
culated by means of diagrammatic technique. Em Analyti- 
cal expressions for negative magnetoresistance have been 
obtained this way for random distribution of scatterer 
in the following cases (i) arbitrary scattering anisotropy 
for low magnetic field B < B tr p where. .B^ = hc/(2el 2 ); 
(ii) isotropic scattering for B 3> B tr .a In the diffusion 
approximation, i. e. when the number of collisions for 
actual trajectories is much greater than unity, this pro- 
cedure givesu 



Ad{b) = Sd(b) - 6d(0) 
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where 7 = t/t V7 b = B/((l + j) 2 B tr ), ip{x) is a digamma 
function, and a is so called prefactor, which is equal to 
unity according to the theory. For x 3> 1 ^(1/2 + %) — 
ln(a;), and the expression (||) coincides with that obtained 
in RpfJ|. The calculations of <5cr(6) beyond the diffusion 
limitoBshow that 5a(b) markedly deviates from this the- 
ory if the number of collisions for actual trajectories is 
not very large. The role of nonbackscattering contribu- 
tion to magnetoconductance has been studied in Ref. [|. 
This contribution has been found to cause the reduction 
of scattering at arbitrary angles and, in contrast to the 
coherent backscattering, the conductivity increasing. In 
the diffusion limit the nonbackscattering contribution is 
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negligible small, but in the case of a strong magnetic field 
B > B tr it should be taken into account. 

It is usual to analyze experimental data by means of 
equation (||) . If this equation describes the magnetic field 
dependence of negative magnetoresistance satisfactorily, 
it is possible to determine t v and its temperature depen- 
dence. 

In our two papers presented back-to-back we put for- 
ward a new approach to calculation and analysis of nega- 
tive magnetoresistance. By representing the quasiproba- 
bility W as a sum of contributions from trajectories with 
given areas we express the negative magnetoresistance in 
terms of area distribution function of closed paths W(S) 
and area dependence of their average lengths L(S). It is 
shown that these are precisely the statistic characteristics 
which can been obtained from the analysis of experimen- 
tal data (see the following paper). In the present paper 
the statistics of closed paths is studied theoretically by 
using computer simulation. This method allows to ob- 
tain the statistic characteristics of closed paths beyond 
the diffusion approximation without any restriction on 
the scattering anisotropy and impurity distribution when 
analytical expressions cannot be derived. 

This paper is organized as follows. In the next section 
we give the necessary formulas and definitions. In Sec- 
tion III the details of simulation procedure are presented. 



The statistics of closed paths obtained from the simula- 
tion is given in Section [Tv| . The results are compared with 
those obtained in the framework of the diffusion theory. 
In Section [y] the magnetic field dependence of negative 
magnetoresistance of the model 2D system are presented 
and analyzed. Both coherent backscattering and non- 
backscattering contributions to magnetoconductance are 
considered. 
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Here, w(S,L)dS gives the density probability of return 
along a trajectory with the length L and area in the in- 
terval (S, S + dS) . Let us introduce the average length L 
of closed paths with a given area in such a way: 



exp 



where 
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L(S) \ 
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W{S)= J —w(S,L). 
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In this case Eq. <M) can be rewritten as follows 



Sa(b) = -2irl 2 G / dS \ W(S) 
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Thus, the area distribution function W(S) and function 
L(S) play a decisive part in the magnetic field depen- 
dence of da. It is clearly seen that these functions can 
be extracted from experimental Sa(b) curves by using 
Fourier transformation. To understand what these statis- 
tic characteristics are let us turn now to theoretical study 
of statistics of closed paths through the computer simu- 
lation of a particle motion over 2D plane with randomly 
distributed scatterers. 



II. BASIC EQUATIONS 

Let us introduce the value wn(S) in such a way that 
ww{S)dS gives the probability density of return after N 
collisions following a trajectory, which encloses the area 
in the range (S, S + dS). In this case Eq. (||) for conduc- 
tivity correction in a magnetic field is written as follows 

Sa(b) = -2nl 2 G a x 

x£;r«.„(W£±3». ( 4) 

n=3 J -°° ^ ' 

In order to take into account inelastic processes destroy- 
ing the phase coherence we include the factor exp(— L/l v ) 
in Eq. (^), where l v is the phase breaking length con- 
nected with t v through the Fermi velocity, l v = UfT v 
and replace the summation over N by integration over 
the path length L. Then, Eq. (||) takes the form 

M&) = -2^ 2 Go£°^{cxp (-^j 



III. SIMULATION DETAILS 

The model 2D system is conceived as a plain with ran- 
domly distributed scattering centers with a given total 
cross-section. It is represented as a lattice M x M. The 
scatterers are placed in a part of the lattice sites with 
the use of a random number generator. We assume that 
a particle moves with a constant velocity along straight 
lines which happen to be terminated by collisions with 
the scatterers. □ The following algorithm has been real- 
ized in computer to obtained the information about the 
statistics of closed paths: 

1. A start point is chosen to coincide with some scat- 
terer near the center of the lattice; 

2. We consider a particle suffers the first collision in 
the start point and begins to move in a random 
direction (see inset in Fig. [I]); 

3. If the particle passes in the vicinity of some scat- 
terer at a distance less than s/2, where s is the 
total cross-section of the scatterer, the collision is 
considered to take place; 
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Then a scattering angle is randomly generated and 
the particle begins to move in the new direction 
until the next collision occurs. The distances be- 
tween two sequential collisions are stored in order 
to calculate the mean free path; 

If the trajectory of the particle passes near the start 
point at the distance less than d/2 (where d is a 
prescribed value, which is small enough), it is per- 
ceived as being closed. Its length L, the number of 
collisions N, the angle 9 between the first and the 
last segments of the trajectory, and the enclosed 
algebraic area, calculated according to 
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(X j+ i - Xj) + 

{x N -Xi), (9) 



where Xj, yj stand for coordinates of j-th collision, 
are kept in memory. Then, the particle continues 
moving in the same direction. Notice that multi- 
returned trajectories are taken into account. Their 
length and area are calculated beginning with the 
start point; 



6. Steps |3j-|5j are repeated until the particle goes out of 
the lattice or a number of collisions exceeds some 
large prescribed value N m ; 

7. Then, another particle is launched from the start 
point in a random direction, and steps are re- 
peated; 

8. After a thousand of starts, a new start point is cho- 
sen, and steps |^-|^ are repeated, for tens of times; 

9. Then a new ensemble of scatterers is randomly gen- 
erated and the whole procedure is done over and 
over again. 

The procedure described allows to find the probability, 
T, for the trajectory to pass at a distance less than d/2 
near the start point, instead of the probability density of 
return, W, which stands in Eq. (Q). But it can be shown 
(see Appendix) that for d <C I these values are connected 
through the simple relationship W = (dl)~ 1 T. 

All the results presented in the paper have been ob- 
tained using the following parameters: the lattice dimen- 
sion is 6500x6500, the total number of scatterers is about 
8 x 10 4 , iV TO = 10 3 , s = 7, and d = 1 (hereafter all the 
lengths and areas are given in units of lattice param- 
eter and lattice parameter squared, respectively). The 
total number of starts, I s , is about 10 6 — 10 7 . The mean 
free path computed for such a system is 77.8. It is close 
to theoretical estimation of mean free path I = (Ngs)^ 1 , 
where N s is a density of scatterers. The calculations were 
carried out also with other parameters M and s differed 
2 — 3 times. It has been shown that this leads only to a 



change in the value of I, but does not influence the depen- 
dence of magnetoconductance on the reduced magnetic 
field b = B/B tr . 

In order to illustrate the 2D system which can corre- 
spond to our model in reality, let us set the lattice con- 
stant equal to 0.5 nm. In this case our model provides 
an example of 2D system with the concentration of scat- 
terers 7.5 x 10 11 cm~ 2 , mean free path I = 38.9 nm, and 
B tr ~ 0.2 T. 
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FIG. 1. The quantity 2nl 2 W N as a function of collision 
number N. Circles are the result of simulation. The line 
is theoretical calculation using Eq. (0). The inset shows 
schematically one of the paths which pass through the start 
point vicinity after N — 8 collisions. 



IV. STATISTICS OF CLOSED PATHS 



A. Simulation results 



It has been shown in Ref. || that if we decompose the 
probability density of return to the origin FT as a sum 
of contributions of paths with TV collisions, Wn> for each 
partial contribution a simple expression is valid: 



2irl 2 W N = {N -2)-\ N > 3. 



(10) 



We emphasize that Ex. ( |To| ) is exact. It has been ob- 
tained for random distribution of scatterers of zeroth ra- 
dius. Corresponding simulation results for our model sys- 
tem are presented in Fig. |IJ As is clearly seen the power 
law works well in the whole range of N. The value of Wn 
is really proportional to (N-2) a with a = (-1.05 ±0.02) 
which is close to the theoretical value a = — 1 . The slight 
difference results from a specific feature of our system. 
The matter is that the scatterers are placed discretely 
rather than continuously as in the theoretical approach. 
They cannot lie closer than one unit cell of the lattice. 
This leads to the fact that unlike the theory a distance 
between two sequential collisions in our model is always 
greater than some critical value of the order of lattice con- 
stant. From this point of view, our model system most 
closely corresponds to real 2D system, where the ionized 
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impurity distribution is discrete and correlated to some 
extend due to Coulomb repulsion at growth temperature. 
It is easy to show in the framework of theoretical ap- 
proach Ref. H that the cutoff of short free paths results 
in more steep Wn(N) relationship as compared with Eq. 
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FIG. 2. Length distribution function of closed paths W(L). 
The solid curve is the result of simulation, the dotted curve 
shows L _1 function. 

The length distribution function for closed paths, 
W(L), is presented in Fig. |[ The function W(L) is de- 
fined in such a way that the value W(L)dL is the proba- 
bility density of return to the origin following a trajectory 
with the length belonging to the interval from L up to 
L + dL. As is seen W(L) is inversely related to trajectory 
length. A drastic decrease of W(L) for L > 7.5 x 10 4 re- 
sults from the restriction on maximal number of collisions 
per one trajectory in our algorithm. As a consequence of 
this restriction the probability of a closed path to have 
the length exceeding the value N m l ~ 8 x 10 4 is negligible 
small. 

Figure ^ shows the area distribution function obtained 
from our numerical simulation. Since the results of cal- 
culations are identical for positive and negative algebraic 
areas, hereafter we present theoretical curves in the pos- 
itive area range only. It is reasonable that W(S) is a 
diminishing function. As is seen it is difficult to describe 
W(S) curve by a power function in the whole range of 
areas. However, for large S the power function S" -1 is a 
good asymptotic of W(S). 

The results of our calculation of L(S) for different val- 
ues of l^p are presented in Fig. f|. As is seen the smaller 
the value of 7, i. e. the greater l v , the greater is the 
average length of closed paths with a given area. For 
large area values S > I 2 , the area dependence of L can 
be described by the power function S 13 with (3 slightly 
depending on the value of 7. As 7 changes from 0.1 to 
0.01, the value of /? varies from 0.55 to 0.62. 




FIG. 3. Area distribution function of closed paths W(S). 
The solid curve is the result of simulation; the dashed and 
dot-dashed curves are the diffusion and improved diffusion ap- 
proximations, respectively; the dotted curve shows 5 1-1 func- 
tion. 



Although the particle motion in our model 2D system 
is a special case of the random walker problem, which is 
studied well enough, to our knowledge there is no ana- 
lytical solution of this problem in a wide range of path 
lengths and areas, involving non-diffusion motion. There 
is no theory which could describe analytically the statis- 
tics of closed paths with small number of collisions. Be- 
low we analyze the simulation results in the framework 
of diffusion theory. 
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FIG. 4. The area dependence of L for different 7 values. 
Solid curves are the simulation results, dashed and dot-dashed 
curves are the results of calculation within the diffusion and 
improved diffusion approximations, respectively (see text). 



B. Comparison with diffusion theory 

Let us now calculate within the diffusion approxima- 
tion the quantities w(S,L), W(S), and L(S) introduced 
in Section O to describe statistical properties of a random 
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walk. A particle trajectory is regarded as being diffusive 
when its length is greater than the mean free path, asso- 
ciated with the transport scattering time. According to 
the diffusion theory the probability density of finding a 
particle in the vicinity of the point r at the moment t, 
P(r,t), does not depend on the initial and final velocity 
directions and obeys the usual diffusion equation 



dP 
~dt 



DAP = 6{t)5(r), 



(11) 



where D is the diffusion coefficient. A solution of the dif- 
fusion equation can be written in the form of the Wiener 
path integral as (see, e.g. Ref. |J) 



P(r,t) 



r(t)=r 



(0)=0 



Vv(t) cxp — 
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(12) 



With help of this formalism we calculate the probability 
density V for a diffusive walk with a constant velocity vf 
and the length L = vpt to enclose the algebraic area S 



V(S, L) 




-L/vf 



dTr 2 (T)0(T) 



P(0,t) 



Here the symbol (— )p(o,t) stands for averaging with the 
Wiener measure dl2) over all closed trajectories. The ex- 
pression for V(S, L) can be obtained by making use of 
relation between this quantity and density of states of a 
fictitious quantum particle with the mass m — h(2D)^ 1 
in a magnetic field:EI 



— ^— cosh 2 
21L 



Tl 



(13) 



When obtaining Eq. (|13|) we have used the expression 
for the diffusion coefficient D — v^t/2, and normalized 
V{S, L) to unity: 



dSV{S,L) = 1. 



(14) 



The function w(S,L), which is connected with V(S,L) 
through the diffusion density probability of return, by 
the following relationship 



<S,L) = 



1 



4nDt 



V{S, L) 



1 



2nlL 



(15) 



is shown in Fig. ^| for different L values. In the same fig- 
ure the results of simulation are presented, too. As would 
be expected, the diffusion theory describes only the 
statistics of long trajectories (see corresponding curves 
labeled 50 1 and 500 I in Fig. ||). For ballistic trajectories, 
there is no agreement between the simulation data and 
the above theory. 

To calculate the function W(S), which is given by the 
Eq. (0), we choose / as the lower cutoff in the integral, 
because the ballistic trajectories are not described within 



the diffusion approximation. Using Eqs. (|T^) and fll5| ) we 
immediately obtain 



ttS 

Anl 2 W{S) = ~ tanh ( — 



(16) 



In Fig. 0, the results of calculations of the function 
4nl 2 W(S) are presented by the dashed line. As is seen 
for large areas S > I 2 , the simple diffusion theory de- 
scribes the simulation results perfectly. For S < I 2 the 
theory gives the values of W smaller than the simulation 
data. 




FIG. 5. The quantity w as a function of S for different L 
values. Solid lines represent the simulation data. Dashed and 
dot-dashed lines show the results of calculation using Eqs. 
( p"3| ) and (|l7|), respectively. Arrows indicate the critical area 
value as it is obtained from the simulation procedure. 

In the framework of the diffusion approximation we 
have calculated the area dependence of the average length 
of trajectories L, introduced by Eq. (^). The results 
of numerical calculation are presented in Fig. |^ by the 
dashed curve. It is evident that the theoretical and sim- 
ulation results are in good agreement only for S > 101 2 . 

Let us suggest some improvement of the above theory, 
which could allow one to describe the statistic of ballistic 
trajectories as well. As is seen from Fig. ^|, at some value 
of S = S c (marked by arrows in Fig. H) each curve ob- 
tained from the simulation procedure reveals a step-like 
behaviour: there are no paths with S > S c . Such a be- 
haviour of w(S, L) is quite clear. A closed path with the 
length L cannot enclose the area larger than the area of 
a circle with the radius L/(2ir). This fact, which has not 
been taken into account above, does not play an essential 
role within the diffusion regime, because the value of Sc 
is much greater than the area enclosed by almost all the 
trajectories. In the ballistic regime, S c is found to be 
equal to the areas enclosed by the most probable trajec- 
tories (see Fig. [| curves corresponding to L = 11, 51). 
Not counting the existence of S c in this case the function 
V(S, L) (and w(S, L) too) is found to be underestimated 
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for S < S c owing to the normalization of V(S,L) ac- 
cording to Eq. (|1J) . The artificial cutoff of the function 
V(S, L) obtained in the framework of diffusion approxi- 
mation (i.e. V{S,L) = for \S\ > L 2 /(4tt) ) and 
following normalization (jl4|) of this function gives 



mi) = ^tanh-(§ 



ttS 



xcosh (tt 9 U?- |6f| 



L 2 



(17) 



where 0(a;) is Heaviside step function. The use of this 
expression allows to describe w(S, L) behaviour for the 
ballistic trajectories much better than Eq. jl^ ) (see the 
dash-dotted curve in Fig. ||). The results of numerical 
calculation of W(S) and L(S) obtained in the framework 
of the improved theory are presented by the dash-dotted 
curves in Figs. || and |j, respectively. A good agreement 
with simulation results is evident for W(S) in the whole 
area range. However, such an approach is too rough to 
describe L-vs-S curves . 

Thus, the theoretical analysis shows that the simula- 
tion procedure works correctly and allows to use it for 
analyzing, among other things, the interference quantum 
corrections to the conductivity. 



V. NEGATIVE MAGNETORESISTANCE 

Before the discussion of negative magnetoresistance let 
us consider the results of calculation of Sa for 6 = 0. This 
will allow us to realize the restrictions of our model in 
view of the finite size of matrix used and limited number 
of collisions per one trajectory and to analyze for what 
7 range the simulation results can be applied to inter- 
pret experimental data for macroscopic samples (when 
the sample dimensions are much greater than the phase 
breaking length). 

To calculate Sa we use Eq. ([!]) assuming that each 
closed path gives a contribution 1 / I s (where I s is a total 
number of paths) to the probability of return. Each con- 
tribution is weighted by the factor exp(— h/lip) to take 
into account the interference distortion by inelastic pro- 
cesses. The final form of the expression for Sa in the case 
of b = looks as follows 



Sa_ 

Go 



2nl 



ex P --T 



(18) 



where summation runs over all closed trajectories. Figure 
^ shows the results of our simulation of Sa as a function 
of 7 = l/ly and, for comparison, the results of theoret- 
ical calculation obtained through the well known exact 
formula 



Sa , . , „ 

_ = ln(l +7 -i). 



As is seen for 7 > 10~ 2 the simulation and theoretical 
data are almost the same. A strong deviation is observed 
for 7 < 10~ 3 where just the trajectories with lengths 
L ~ L v = 7 _1 1 > N m I have to give an essential contri- 
bution to Sa value. These trajectories are not considered 
in our model. Thus, we believe that for 7 = l/l v > 10~ 2 
our model system is equivalent to an unbounded 2D sys- 
tem and the simulation gives correct results. 
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FIG. 6. The interference quantum correction to the con- 
ductivity in zero magnetic field as a function of 7. The solid 
curve shows the simulation results, the dashed line is the re- 
sults of calculation according Eq. ([l9|). 



A. Backscattering contribution 

Now we are in position to discuss the magnetoconduc- 
tance anomaly due to suppression of quantum interfer- 
ence corrections by a magnetic field. To calculate a mag- 
netic field dependence of Sa for our model system we 
follow the ordinary way: the contribution of each closed 
trajectory to Sa is multiplied by the factor, which allows 
for the interference distortion by the magnetic field. The 
final expression for Sa(b) takes the form: 



8a{b) 
Go 



2tt/ x - 



cos 



(l+7) 2 ^ 



exp 



(19) 



This formula is not meaningful from the experimental 
point of view, because the absolute value of quantum 
correction is not a measurable quantity. The magneto- 
conductance Aa(b) = a(b) — a(0) = 5a(b) — Sa(0) is 
usually experimentally available. 

The results of calculation of Aa(b) for the model sys- 
tem are presented in Fig. [?] for several 7 values. They are 
in excellent agreement with the results of numerical cal- 
culations carried out beyond the diffusion approximation 
in Ref. ||. This lends support to the correctness of pa- 
rameters chosen for our model system: the set of param- 
eters used turns out to be suitable to simulate adequately 
Aa(b) in the ranges of b and 7 under consideration. 
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FIG. 7. The magnetic field dependence of Act for different 
7 values. Solid curves are the results of simulation, dashed 
curves are the diffusion limit ([]). Only backscattering contri- 
bution is taken into account. 




FIG. 8. The magnetic field dependence of A a for 7 = 0.1. 
Solid curves are the results of simulation. Dashed curves are 
the results of fitting procedure with the use of Eq. ^ made 
within two different ranges of magnetic field: b < 1 (a) and 
b < 5 (b). The fitting parameters are the following: a = 0.49, 
7/ = 0.087 (a) and a = 0.65, 7/ = 0.12 (b). 

In the same figure the results of calculation of Aa(b) 
in the framework of diffusion approximation with the ex- 
pression (H) are presented too. It is seen that this formula 
does not describe the simulated Aa(b) curves. Even in 
the case of 7 = 0.01, when the condition 7 <C 1 is seem- 
ingly fulfilled and the diffusion theory should work well, 
the formula (||) gives the value of Act substantially larger 
than the simulation data: for b = 0.1 the difference is 
about 25%. 

As discussed above (see Section |) Eq. (||) is widely 
used by experimenters to extract the values of 7 from 
experimental data. Two fitting parameters are available 
in such a data processing: just as the value of 7 so the 



prefactor a. To check the validity of this method in the 
case when the diffusion approximation does not work, we 
have performed the fitting procedure for Act(6) curves 
simulated. In Fig. || the results of fitting procedure of 
Act (b) made for 7 = 0.1 within two different ranges of 
magnetic field b < 1, and b < 5 arc presented by dotted 
and dashed curves, respectively. As is seen the simulated 
curve is best described by Eq. (||) in low magnetic field 
range. However in both cases the ratio t/t v found from 
the fitting procedure and designated as 7/ slightly differs 
from the value of 7 used in the simulation. The differ- 
ence is about 10% for b < 1 fitting range and 20% for 
b < 5. The value of prefactor is found to be significantly 
less than unity. 

The results of such a data treatment for several 7 val- 
ues are presented in Fig. |^. It is clearly seen that the 
difference 10.. 30% between 7/ and 7 takes place for all 7 
values. The value of prefactor a is always less than unity, 
and increases with decreasing 7. 
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Y Y 
FIG. 9. The parameters a (a) and 7/ (b) as a function 
of 7 value fed into the simulation. Open and solid cir- 
cles are the results of fitting of Aer(&Y, obtained without 
non-backscattering contribution, to Eq. (g) in the ranges b < 1 
and b < 5, respectively. Solid lines are the guide for an eye, 
the dashed line shows 7 = 7/ function. 



B. Nonbackscattering contribution 

Up to this point we considered the coherent backscat- 
tering correction to conductivity. The coherent paths for 
backscattering contribution are schematically depicted in 
Fig. |To|a. In this case the interference takes place when 
the points A, 1, N and B are close to one line. However 
as it is shown in Ref . ^ there is one more variant of coher- 
ent paths for the same configuration of scatterers, when 
the finish point B lies close to the line passing through 
the points 1 and 2 (see Fig. |To|b). In the latter case the 
one of the two interfering waves is scattered twice by the 
scatterer 1 (which is regarded in this paper as a start- 
ing point), while the other one passes the point 1 with- 
out scattering. This provides so called nonbackscattering 



7 



contribution to conductivity. It was shown inQ that both 
contributions are related to the probability of return to 
the starting point for a given trajectory, the conductivity 
correction due to nonbackscattering processes is positive. 
In our case such processes can be easily taken into ac- 
count by multiplying each term in Eq. (20) by the factor 
(1 — cos(0j)), where 0i is the angle between the first and 
last segments of the i-th closed path: 



Ao-(fr) 



2ttZ 
'Zd 



]T[...](l-cos(^)). 



(21) 



simulated data as experimental ones, i.e. Aa(b) curves 
have been fitted to the calculated with Eq. (||) values. 
The fitting results are presented in Fig. O. As one would 
expect, the values of a and 7/ obtained from the low mag- 
netic field fitting (b < 1) are very close to corresponding 
values in Fig. |9|. This obviously results from the fact 
that the solid and dashed curves in Fig. [ll] are close to- 
gether for these b values. As for the parameters obtained 
from the whole magnetic field range fitting, taking into 
account of the nonbackscattering contribution results in 
decreasing both a and 7/ values. 





FIG. 10. Two types of coherent paths and corresponding 
diagrams relevant in the first order in (kpl)~ , which are 
responsible for backscattering (a) and nonbackscattering (b) 
contributions to weak localization, for given configuration of 
scatterers. 

Shown in Fig. [H] are the magnetic field dependencies 
of Act, calculated for different 7 values with the help of 
Eq. fl21j). For comparison, the results computed with the 
formula (^0|) are presented too. It is clearly seen that 
both types of calculations give close results in the low 
magnetic field range b < 0.5. For higher magnetic fields 
the inclusion of nonbackscattering contribution leads to 
a decrease in magnetoconductance. Such a behaviour of 
the magnetoconductance is in a good agreement with the 
results of numerical calculations presented in Ref. ||. 

As in the previous subsection we have attacked the 



<! 




FIG. 11. The magnetic field dependence of Act for differ- 
ent 7 values. Dashed and solid curves are the results of cal- 
culation without and with non-backscattering contribution, 
respectively. 




10" 10 10 10 

Y Y 
FIG. 12. The parameters a (a) and 7/ (b) as a function of 
7 value fed into the simulation. Non-backscattering contribu- 
tion is taken into account. All designations are the same as 
in Fig. §. 



Thus, the standard method of experimental data pro- 
cessing with the use of Eq. (^) allows to extract the phase 
breaking time (or length) with the accuracy (10 — 30) % 
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in wide ranges of 7 and b. The prefactor a obtained this 
way is found to take the values from 0.3 to 0.7, i.e. sub- 
stantially less than unity. This reveals a possible reason 
of small value of prefactor obtained from different exper- 
imental data treatment. This need not be a consequence 
of an electron-electron interaction as discussed in some 
papers, but can result from the fact that in real 2D sys- 
tems the rigorous conditions of diffusion approximation 
are not fulfilled even in the case of rather small 7 and b. 



VI. CONCLUSION 



This paper is intended to demonstrate the possibility of 
obtaining the information about the statistics of closed 
paths from the analysis of anomalous magnetoconduc- 
tance in 2D systems. In particular we have shown ex- 
plicitly that the statistic characteristics such as the area 
distribution function of closed paths and area dependence 
of their average lengths determine the magnetic field de- 
pendence of magnetoconductance. These functions have 
been studied by using the computer simulation method. 
It has been shown that in the ballistic regime the area 
distribution function of closed paths deviates from S 1-1 
power law, which holds in the diffusion regime. The theo- 
retical analysis of simulation results has shown that such 
a behaviour of W(S)l is mainly connected with the ex- 
istence of critical area value, which corresponds to the 
maximal area enclosed by a path with a fixed length. It 
has been shown that the area dependence of the average 
length of closed paths, introduced through Eq. (|^), can 
be well described by the power function, L cx S 13 , with (3 
varying in the range 0.55 — 0.62, when 7 = r/r^ changes 
from 0.1 to 0.01. 

The results of simulation have been used to calcu- 
late the magnetoresistance of the model 2D system. 
Both backscattering and nonbackscattering processes 
have been taken into account. The calculated magnetic 
field dependencies of Act have been processed as experi- 
mental ones by a standard manner with the help of Eq. 
(||). It has been shown that the fitting procedure gives 
t /t v ratio, which differs from that used in the simulation 
by the factor 0.8 — 1.3. The value of prefactor obtained 
this way is always less than unity. 
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APPENDIX 

For the sake of simplicity we suppose here that the 
start point coincides with the origin. For the value of 
probability T that a trajectory passes near the origin at 
the distance less than some prescribed value d/2 we can 
write 



T=^T Ar , 

N=3 



(22) 



where T/v is the probability of return after N collisions. 
The value of T/v can be expressed through the probabil- 
ity density to experience the (N — l)-th collision in the 
point r, Wn-iM; 



T N = J drV{r)W N -i{r), 



(23) 



where V(r) is the probability that the particle moves 
without collisions from the point r to the circle of the 
radius d/2 around the origin. It is easy to show that 



x 1 ( d \ ( r 

Vyr) — — arctan — exp — - 
7r \ 2r I \ I 



(24) 



Then, for T/v we have 



rr- 1 f , ( d 

T/v = — I dr arctan — I exp 
7T ./ V 2r 



(-y)wV-i(r). (25) 



For d <^ I, the inverse tangent in Eq. (|25]) can be replaced 
by its argument. Analysis shows that for d/l = 10~ 2 (it 
is true for our case) this gives an error in calculation of 
Tjv less than 0.1 % for N > 3. Thus the expression for 
the probability T becomes 



N=3 J 



(26) 



Comparing Eq. Jm&) with the expression for the proba- 
bility density WaW 

1 Nm r d 

w =^EjT acp (ri) w ^ (27) 

AT — Q J 



JV=3 

we can conclude that 



W ~ (Id^T. 



(28) 
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